Two-dimensional quantum spin-1/2 Heisenberg model with competing interactions 
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We study the quantum spin-1/2 Heisenberg model in two dimensions, interacting through a 
nearest-neighbor antiferromagnetic exchange (J) and a ferromagnetic dipolar-like interaction (Jd), 
using double-time Green's function, decoupled within the random phase approximation (RPA). 
We obtain the dependence of ksT^jJd as a function of frustration parameter 8, where T c is the 
ferromagnetic (F) transition temperature and 5 is the ratio between the strengths of the exchange 
and dipolar interaction (i.e., 8 = J/ Jd)- The transition temperature between the F and paramagnetic 
phases decreases with S, as expected, but goes to zero at a finite value of this parameter, namely 
8 — S c = 7r/8. At T = (quantum phase transition), we analyze the critical parameter S c (p) 
for the general case of an exchange interaction in the form Jij = Jd/r^, where ferromagnetic and 
antiferromagnetic phases are present. 



I. INTRODUCTION 

Considerable attention during the last decade has been 
devoted to the investigation of systems with long-range 
interactions. In particular, the interest in films and quasi- 
two-dimensional systems have attracted attention mainly 
due to their technological applications as, for example, in 
electronics, data storage, catalysis in the case of metal- 
on-metal films, biotechnology, and pharmacology in the 
case of molecular films. The magnetic size unit and 
its thermal stability are essential points to be consid- 
ered in order to obtain a good performance in magneto- 
optical recording^]. The magnetic properties of these 
films depend on the subtle interplay between the long- 
range antiferromagnetic dipolar interaction, the short- 
range rotational invariant exchange, and the magnetic 
surface anisotropy. The presence of antiferromagnetic 
domains (AF) was observed in epitaxial thin films by 
the polarization-dependent X-ray magnetic linear dichro- 
ism spectra microscopy |2|. A possible intrinsic mecha- 
nism leading to AF domains is the competition between 
anisotropy and dipolar interaction 

In two-dimensional lattices and for spins with rota- 
tional symmetry, long-range order does not occur at any 
finite temperature, for quantum and classical models 
with short-range interactions^. That is the case when 
only short-range exchange interactions are present, for 
the XY and Heisenberg models, for example. For all real 
systems there is a long-range dipolar interaction, with the 
important property that it breaks the symmetry between 
out-of-plane orientation of the spins and the in-plane ori- 
entation of the spins in the ordered state. The presence 
of dipolar interactions in systems with rotational sym- 
metry (exchange) may stabilize long-range order at finite 
temperatures, for both classical and quantum models[5j. 

The competition between long-range antiferromagnetic 
dipolar and short-range ferromagnetic exchange inter- 
actions in two-dimensional uniaxial (Ising) spin sys- 
tems is responsible for a very rich phcnomenologi- 
cal scenario concerning both their equilibrium statisti- 



cal mechanics l a LZt ISt l a and non-equilibrium dynamic 
propertiesfll lirllll PlTl Hf| . By means of Monte Carlo 
simulations and analytical calculations of the ground 
state, Maclsaac et al.0 studied the two-dimensional 
spin-1/2 Ising model with ferromagnetic exchange and 
antiferromagnetic dipolar interactions, and have shown 
that for S = J j Jd < 0.425, J and Jd being the strength of 
exchange and dipolar interactions, respectively, the anti- 
ferromagnetic state is stable. For 5 > 0.425 the AF state 
becomes unstable with respect to the formation of striped 
domain structures, i.e., to state configurations with spins 
aligned along a particular axis forming a ferromagnetic 
stripe of constant width h, forming a superlattice in the 
direction perpendicular to the stripes. Monte Carlo re- 
sults at low temperatures give further support to this 
proposal, at least for intermediate values of 5. 

There are few numerical results concerning the equi- 
librium statistical mechanics of quantum models with 
competitive exchange and dipolar interactions in two- 
dimensional lattices [l4lll5| . The two dimensional antifer- 
romagnetic (AF) Heisenberg model has been investigated 
by many authors to explain, for example, the magnetic 
mechanism of high T c superconductivity |l fl Anti- 
ferromagnetic fluctuations are believed to play an im- 
portant role in the superconductivity of the cuprates|l7|. 
such as La2Cu04, which is well described by a quasi- two- 
dimensional quantum spin-1/2 Heisenberg AF model. 
Chandra and Doucot|l£| have studied the square-lattice 
Heisenberg model at T = with next-nearest-exchange 
coupling and suggested that the AF order is destroyed 
due to the competition between the nearest and next- 
nearest exchange interactions. Although frustration ef- 
fects are very effective to destroy the AF order in general, 
they may not be effective in high T c materials, consider- 
ing the extended nature of the holes which destroy the 
AF order. A good example of quasi- two-dimensional AF 
materials in which the dipolar interaction (with presence 
of frustration) is comparable to the exchange coupling 
are the so called high-T c superconductors RBa2Cu307_ a; 
(where R stands for rare earth) 0] . 
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^From a theoretical point of view, the double-time 
Green's function (GF) theory[2(| is both a convenient 
and an effective theoretical framework for interpretation 
and forecasting of various characteristics of matter at 
all temperatures. The development of approximation 
schemes for the GF approach has focused on decoupling 
their equations of motion. This decoupling is usually 
chosen for convenience or for reasons which are essen- 
tially ad hoc. The consistency of the basic decoupling 
approximation with relevant operator identities is not al- 
ways assured. Several first-order decoupling have been 
proposed in the literature: firstly, the random phase ap- 
proximation (RPA) was applied to the quantum spin- 
1/2 Heisenberg ferromagnetic [l^, and extended to in- 
clude a long-range interaction J(r), which depends on 
the distance r between spins as an inverse power law 
J(r) = J/r p , where p > d (dimension) HJ |H . This 
RPA decoupling predicts the absence of magnetic order 
in finite temperature and low-dimensionality (d = 1,2) 
for p > 2d, in accordance with the generalized Mermin 
and Wagner theorem[24|]. The quantitative results of the 
phase diagram in the (T — p) plane in the d < p < 2d 
region obtained by RPA |2^| are in accordance with quan- 
tum Monte Carlo simulations [1^. 

The aim of this work is to investigate the results of the 
competition between the exchange and dipolar interac- 
tions in the two-dimensional quantum spin-1/2 Heisen- 
berg model. In Sec. II, the model is presented and 
treated by the double-time GF technique in the RPA 
decoupling. In Sec. Ill, we analyze the phase diagram in 
the T — S plane. Conclusions and remarks are summa- 
rized in Sec. IV. 



II. MODEL AND GREEN'S FUNCTION 
METHOD 



ferromagnetic order. Consequently the ground state of a 
system determined by the dipolar interaction alone dif- 
fers from the ground state of a system determined by the 
exchange interaction alone. In the absence of the dipo- 
lar interaction (i.e., J a — 0), the isotropic Heisenberg 
model is regained and in two-dimensional lattices it does 
not present long-range order at T > (i.e., T c = 0). 
When both interactions are present the system is inher- 
ently frustrated. 

The double-time Green's function {(A(t); B(0))) is de- 
fined by[U 



((A(ty,B(0))} = -iO(t)([A(t),B(0)]), 



(2) 



where 9{t) is the step function, [A, B] is the commutator 
of operators A and B, and (■ ■ •■) denotes an average with 
respect to the canonical density matrix of the system at 
temperature T. The time-Fourier transform of Eq. (2), 
{(A; B)) E , satisfies the following equation of motion 



E((A;B)) E 
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([AB]) + (([A,H];B)) E , (3) 



and ({[A, 7i] ; B)) E obeys an equation similar to Eq. (3), 
with a higher-order Green function appearing on the 
right side. In this way, an infinite set of coupled equations 
is generated. Therefore, an approximation (decoupling) 
is used to obtain the Green function. 

The correlation function (BA(t)) is obtained by the 
spectral representation theory, which gives 



(BA(t)} = I J{w)e- lwt dw = 

i {G(w + ie) - G(w - ie)}e- iwt dw 
eP w - 1 



(4) 



In order to study the consequences of both quantum 
effects and frustration, we propose here a spin model to 
describe the destruction of the F order, represented by 
the following Hamiltonian 



H = J J2 S ' 

<i,3> 



\ - Si ■ Sj 

(id) 13 



(1) 



where S, = (Sf,Sf,Sf) is the spin-1/2 operator at site 
i. The first term, J, is the strength of the exchange 
interaction, and the sum is over all nearest-neighbor 
pairs, < i,j >, on a square lattice. The second term, 
Jd = (g[i B ) 2 /a 3 (g is the Lande factor, [i B the Bohr mag- 
neton, a the lattice constant), represents a long-ranged 
dipole-dipole interaction and the sum is over all possible 
pairs of atoms on the square lattice. The dipolar interac- 
tion tends to align the spins in the F system at low tem- 
peratures (T < T c ), while the exchange interaction (an- 
tiferromagnetic, J > 0) tends to destroy the long-range 



where G{E) = ((A;B)) E and s -> 0. Note that Eq. 
(4) is the required spectral representation for the time 
correlation function, where J(w) is the spectral intensity 
of the function (BA(t)) (in fact, its Fourier transform 

H). 

With A = 5+ and B = S l , where the spin operators 
are defined by the usual commutation rules, we obtain 
fromEq. (3): 



E((Sp,S r )) E = -5 gl +2j2JA((S;S+;S r )) E 

-((S*Sp,Sr)) E> (5) 

where m = \Sg) is the magnetization per spin, Jji = — J 
for nearest-neighbor sites and Jji = 4r for the dipolar 

interaction (with S = J/ ' Jd). 

The key problem of a first-order decoupling pro- 
cedure is essentially to express the Green's function 
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((S^S^; Sf)) _ in terms of lower-order Green's func- 
tions, which enables one to solve the infinite chain of 
equations of motion in an approximate way, and that 
can be expressed in the following form: 

((S z a S+;Sf)) E cm((S+;Sf)) E . (6) 

Using the decoupling (6) in Eq. (5), we obtain 



EG gl (E) = -S gl + 2mJ2Jj a {G 9 i(E) - G jt (E)} (7) 



where G gl (E) = ((S+;Sf)) E . 

The method of calculation we use is the same as 
Nakano and Takahashi[23j. and we shall not reproduce 
the details. Following their notation, we obtain from 
Eq. (7) the Fourier transform of the Green's function 
G K {E) = F{G g i{E)}, defined by 



G K (E) = J2G g i(E)e-^ 

9,1 



and, in this way, we find 
G K (E) = 



tt(E — Ek) 
where the magnon energy Ek is 

E k = 2m( J - J & ), 

with 

^ = E^ e " k,(r9 " r!) - 

(9,1) 



(8) 

(9) 
(10) 
(11) 



Using the low-k expression in two dimensions (as in 
Ref. 23) for p — 3 (the ferromagnetic interaction is as- 
sumed to decay with the distance between spins, r, as 
J(r) = J/r p ) we obtain: 



E k = 2mJ d k(ir 2 /8 - 6k). 



(12) 



Using Eqs. (4), (8) and (9), we obtain the correlation 
function (S-S+); for spin S = 1/2 we have (S~S+) = 
1/2 — m and, therefore, the magnetization is written in 
the form: 



l]Tcoth(/3£y2) 



(13) 



III. RESULTS AND DISCUSSION 

In the limit m — > 0, we find, from Eqs. (12) and (13), 
the critical temperature (T c ), which is given by: 



k B T c 
Jd 



i 



N ^ k(ir 2 /8-6k) 



(14) 



FIG. 1: Dependence of the reduced critical temperature 
k,BT c / Jd on the parameter S for the quantum spin-1/2 Heisen- 
berg model with ferromagnetic and dipolar interactions with 
p = 3. The insert is the behavior of 5 c (p) as a function of the 
parameter p 



In the thermodynamic limit (N — > oo) one must replace 



the sum jj ^(k) by an integral 



(27r) d 



iP(k)d d k in 



1BZ 



the d-dimensional k-space, where 1BZ denotes the first 
Brillouin zone. Then, the integral (14) is obtained by 
using the expansion (12), and the critical temperature is 
given by 



k B T c 



4ttS 



In 



(15) 



Note that the critical temperature vanishes at 5 C , which, 
within the present approximation, is given by S c = tt/8 ~ 
0.39. 

Numerical results for T c are shown in Fig. 1. The 
parameter S is a measure of the strengthen of the frus- 
tration; therefore, long-range order decreases as <5 in- 
creases. Consequently, T c {8) goes to zero at the critical 
value S = S c — tt/8. When the interactions have the 
general form Jd/r p , the existence of a finite-temperature 
phase transition occurs for 2 < p < 4 in two dimen- 
sions. By numerically performing the sum in Eq. (13) 
for 2 < p < 4 (see Ref. 23 for the general dispersion re- 
lation), we obtain the critical parameter S c (p) as a func- 
tion of p at T — (see the insert of Figure 1), where 
above the curve the ordered phase is antiferromagnetic, 
while below it the long-range order is ferromagnetic. We 
note that 5 c (p) increases as p decreases, and when p ap- 
proaches p = 2 we have a divergence (T c also diverge) 
in S c (p = 2) and around p ~ 3.6 a minimum point ap- 
pears. The point p = 4 corresponds the to Kosterlitz- 
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Thouless phase transition |23j, and with the presence of 
the AF nearest-neighbor interaction the critical temper- 
ature vanishes at 8 c (p = 4) ~ 0.793. For p > 4, the 
critical frustration parameter, 5 c (j>), is zero, as physically 
expected. 

IV. CONCLUSIONS 

We study the phase diagram of the quantum spin-1/2 
Heisenberg model with competing interactions (i.e., 
presence of AF exchange and dipolar interactions). The 
influence of the frustration is analyzed through the 
variation of the parameter <5 = J/Jd- We observe that 
T c {5) is null when we reach the critical value S c = tt/8. 
We also analyzed the case of the general dependence of 
the long-range interaction in the form Jjj = Jd/f^j, arid 
verified that T c tends to zero at the critical parameter 
S c (p). The dependence of S c (p) on p (2 < p < 4) present 
a divergence at p = 2 (p < 2 being the non-extensive 
regime) and a minimum point around p ~ 3.6, with 
finite value at p — 4 (Kosterlitz- Thouless transition). 
Although quantitative results are not expected to be 
very precise, within the approximation we applied, we 
believe that we obtain reliable qualitative results for the 
whole range of the parameters of the Hamiltonian. In 
particular, our results agree with previous ones, when 
available. 
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